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MCTM - C0 2 comprises 95% of the composition of the “f f'rom leTthan 

the martian atmosphere also has a high wi t h the occurrence of global dust 

0.2 to greater than 3.0, primarily on a seasonal b « sls 0 tQ greater than 1. with 

storms during southern spring ^ [2] _ lH ,%V C0 2 is an active absorber and emitter 


Storms <iUXXU£ w^..— 0 i , _ . , r o 1 _ 

large amounts occurring at winter polar ®f l J * ^2 x " provide significant heating 
in near-IR and IR wavelengths; the near-IR Howevef. dust and ice- 
of the atmosphere, and the 15 band prov P wavelengths, and are highly absorbing at 
ol.ua aerosols have high ...Ctoylng albedo.. ta ~ " LTTStel. «n.rgy budg.C [5). 


cloud aerosol, have high soatterlng .lo.ooe. in .^herlc energy budget IH 
infrared wavelengths, and are as importan 2 simultaneously in a model is difficult. 

Including both C0 2 and aerosol ra a ^e scattering code while C0 2 radiative transfer 
:osol radiative transfer requires a muitlple- scattering • The f eblelB can be so lved 

ictly & by ^^ertlng 6 * hi^^^^absoi^tMc^for ^each spectral line into a 

t.dlJ“ t ,nd h '.^ir d ST.^u”dlSUTdr.^S aerosol -».U 


ceaiuub emu. ~ • — — _ 

scattering and the banded absorption structure or 
C0o that prompts most radiative- transfer studies 
of the martian atmosphere to consider either a 
pure-C0n or pure-dust atmosphere. This approxi- 
mation simplifies treatment, but is inaccurate. 

One alternative technique that has recently 
been developed for atmospheric applications is 
the exponential -sum or k-distribution approxima- 
tion [8-21]. The transmission of a homogeneous 
atmosphere is actually independent of the order- 
ing of the absorption coefficient, k, in frequen- 
cy space within a spectral interval, depending 
only upon the percentage of the spectral interval 
that has a particular value of k. The percentage 
of the spectral interval which has values between 
k and k + Ak can be formulated in a probability 
density function f(k) shown schematically in Fig. 
2 The chief advantage of the exponential -sum 
approach is that the integration over k space of 
f(k) can be computed more quickly than the inte- 
gration of k^ over frequency. The exponential- 
sum approach is superior to the photon- 
path- distribution and emissivity tech- 
niques for dusty conditions [22, 19, 

23] . Our work is the first application 
of the exponential -sum approach to mar- 
tian conditions. 

THEORETICAL APPROACH . The transmit- 
tance of the 15 jm band and the near-IR 
bands of C0 2 was computed using the 
FASC0D line-by-line transmittance model 
[24] , modified for martian condi- 
tions. Computations with the modified 
FASC0D model were made at 3 tempera- 
tures (125K, 200K, 300K) and 5 pres- 
sures (100 mb, 10 mb, 1 mb, 0.1 mb, 

0.01 mb); these cover the range of tem- 
perature and pressure currently observ- 
ed in the atmosphere at all latitudes, 
seasons, and altitudes up to 40 km, and 
also can be used for early Mars, dense 
atmosphere studies. The near-IR and 
15/an bands were broken into spectral 
sub-intervals (see Fig. 1), and the 
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temperatures of 200 and 300K [25] . 
Also shown are the sub- intervals used 
for the 8 term and 16 term fits. 



Figure 2. A schematic illustration portraying 
the essence of the exponential -stun approach. 

(a) shows a schematic of absorption line spectra 
at two different pressures. In (b) the two pro- 
bability density functions f(k) associated with 
(a) are illustrated. The shaded area depicts 
the strongest absorption (i.e., largest k) for 
the same spectral interval (i.e., f(k) for dif- 
ferent pressures are correlated) . Integration 
of f (k) over k replaces the integration of k^ 
over v (modified from [28], [29]). 
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( 1 ) 


transmittance for each layer, Tr, as a function of C0 2 column abundance within that layer, u, 
was fit with a series of weighting coefficients, a^, and exponential coefficients, b^. 
n 

Tr(u) « 2 a^exp( -b^u) 

We tried fitting 1 procedures based on Wiscombe and Evans [14] and an improved version of Acker- 
man et al. [13], and found both yielded similar results. Both of these procedures avoid the 
ill-conditioning of earlier exponential- sum routines, and produce more accurate and unique so- 
lutions [13,14], The exponential -sum fit reproduced the FASCOD transmittance s to better than 
10'* for all CO, abundances considered. The frequency sub-intervals are picked to try to 
minimize the variation in line strengths within the sub-interval. As shown in Fig. 1, 2 sub- 
intervals covered the band center, 2 covered the far line wings, 2 covered the near line 
wings, and 2 covered the transition from wings to band center. The vertical inhomogeneity o 
the atmosphere is treated by using homogeneous layers, with the absorption coefficients or 
all layers correlated in frequency space, i.e. the sub - intervals of the spectral b ^ nd '^ h 
have the maximum absorption also have the largest k values (Ackerman et al. . 1976) An inter- 
polation is used for temperatures and pressures which fall in between the values at ^ich the 
exponential -sum coefficients were computed. We compared a logarithmic interpolation 
linear interpolation, and found the logarithmic interpolation to be more accurate 

INCORPORATION IN Min .TIPLE- SCATTERING MODELS . Vertical optical depths of C0 2 absorption tor 
each term number, frequency sub- interval , and atmospheric layer are bjU. C0 2 is combined with 
dust and cloud in that the total optical depth T t , single -scattering albedo and phase 

function P for each term, frequency su^- interval , and layer are given by[16]: 

T. - b.u + r + % + r a + r s + r a 
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(4) 


where: r R - Rayleigh scattering optical depth for that layer 

r D » C - Dust (D) and Cloud (C) scattering (s) and absorption (a) optical depth for 

’ that layer 

P D » C > R - dust (D), cloud (C), and Rayleigh scattering (R) phase function 
The multiple - scattering code is run once for each term in the sum using the T t w i$ and P 
appropriate to that term, and the resultant fluxes, F,. (or intensities) are then summed and 
weighted by ^ to give the total flux, F, (or intensity! over the frequency sub-interval: 


F - 2 a. F. (b.u) (5) 

i-1 1 1 1 

Tables of the exponential -sum coefficients, a^ and b^, 
can be obtained from the authors . 

NUMERICAL STUDIES . The number of terms in the 
series of exponentials (n in equation 1) and the num- 
ber of sub - intervals into which the spectral band is 
broken can be varied to increase the desired accur- 
acy. Figure 3 shows a comparison between the 15 pm 
cooling rates computed using various numbers of fre- 
quency sub- intervals and terms within each sub- inter- 
val. The temperature of the atmospheric model rose 
linearly from 150K at the surface to 160K at 10km 
altitude, and then fell linearly to 130K at 40km alti- 
tude. These temperatures are typical of the winter 
polar atmosphere on Mars. Using 4 terms and 4 sub- 
intervals results in an error of less than 10% in the 
lowest 10 km, but results in substantial errors at 
higher altitudes. The explanation is that the 4 term 
fit could not capture the effects of both the band 
center and the wings at all altitudes, and in this 
case we emphasized the wings, which are more important 
at lower altitudes where the band center is saturated 
(see Fig. 1). Adding additional terms from 8 to 16 is 
far less noticeable. 



Figure 3. 15 ^m band cooling rates 
calculated with an exponential -sum 
approach using 4 frequency sub- in- 
tervals and 4 terms in the sum, 10 
sub- intervals and 8 terms, and 10 
sub- intervals and 16 terms. A win- 
ter polar temperature profile is 
used (see text), but only 200K expo- 
nential-sum coefficients are used 
for comparison purposes. 
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We also separately derived exponential -sum coeffi- 
cients based on the Gal' tsev and Osipov [25] line-by- 
line calculations of the 15 /im band. The number of 
terms in the fit (n in equation 1) were varied, to 
confirm our FASCOD work on the number of terms requir- 
ed to give an accurate fit. Figure 4 shows cooling 
rates computed with several fits to the Gal'tsev and 
Osipov parameterizations , as well as to the FASCOD 
transmit tances . Cooling rates compare favorably in 
the lower atmosphere (below 10 km altitude). Gal'tsev 
and Osipov only considered temperatures as cold as 
200K. We extrapolated the temperature dependence to 
colder temperatures [26]. Discrepancies between the 
Gal tsev and Osipov and FASCOD cooling rates are due 
to the inaccuracy in the temperature extrapolation, 
particularly at the higher, colder altitudes. Again, 
note that the 4 and 5 term fits become inaccurate 



above 10 km altitude, as in Fig. 3. Also, the 15 term 
fit yields no marked improvement over the 8 term fit, 
as also shown with the FASCOD fits in Fig. 3. 

Cooling rates computed using the FASCOD exponen- 
tial-sum transmit tances also compared well below 10 km 
altitude with cooling rates computed using the Pollack 
et al. [27,4] parameterizations of C0 2 transmit- 
tance. At higher altitudes, discrepancies exist due 
to the use of the strong- line approximation by Pollack 
et al, which emphasizes the effects of the line 
wings. However, comparisons of our exponential -sum 
transmittances with techniques commonly used for the 
terrestrial atmosphere have indicated that our exponential - sum transmittances could be in 
error above 20 km altitude. J. Pollack is seeking to modify our exponential -sum approach to 
use two sets of sums, one applied to the line centers and one applied to the line wings. 
Initial tests seem to show better agreement in the upper atmosphere of Mars. 


Figure 4. 15 ^m band cooling rates 

calculated with 4 types of exponen- 
tial-sum coefficients: 4 and 15 

term fits to the transmission para- 
meterizations of [25] (abbrev. G&O) , 
and 4 and 8 term fits to the FASCOD 
transmission model. The winter po- 
lar atmosphere model is used, with 
logarithmic temperature interpola- 
tion for the exponential coeffi- 
cients . 


B.L. Lindner acknowledges support by NASA contract NASW-4444. 
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MAGMATIC PROCESSES ON MARS: INSIGHTS FROM SNC METEORITES; J. 
Longhi, Lamont-Doherty Geological Observatory, Palisades, NY 10964 . 

The composition and petrology of the SNC (shergottites-nakhlites-Chassigny) meteorites reveal a 

g?? S2 l m i» 

f . Vn crV* p- low A1 liauids (2) As might be anticipated, calculated densities of these liquids are high (2.75 

MORB's end comment AotaL m which plegioclese ctystellte eterly. These low-Al .S'?, 

S Tie for the Acmes. Thts depletion m.y be eceomplished by the .ccumql.non of olt.tne end pyroxene 

“tL IT£2£ of SNC trace element and isotopic composition. Figute 1A 

shows the REE concentrations of two calculated parent liquid compositions for Nakhla. Chasstgny parent ltqut 
Srstaflar pL™ S“. « unused rock costing of large cumulus augite ami minor oh.me c ? s.ds set m 
* n;_. motriv The ’closed system' calculation assumes that no net changes in the intercumu 
Hqu^Jok^lace after accumulation of the pyroxene, but does allow for partial equilibration of the cumujus ciystals 
nnH framed liouid The open system' calculation is a direct calculation based upon the composition of Nakhla augn 
^ Sfi Zlkion coeffi^ents of (7); this calculation allows for the possibility of migration of mtercumulus 
liquid. B^th patterns are similar and show dramatic Ught REE enrichment. This fractionation is truly remarkable m 

light of the ENd v ^e of +16 (6) which requires that the source had a long term pattern of light REE de P e 10 ’ , 
s u- • -7 ^ upT a 7 orwn A nattem in Fig 2B Compounding the situation is the low A1 content of the 

^r^xs ss ^”r^.T(™ i rm F .‘ £ Ai^ e Smpo„L,. r a, 4“^!“ 

garnet which is the most effective REE fractionating agent, cannot have been a residual P h “* in **£“*"} 
magma’s source region; niether is there much allowance for removal of augite at low pressure. The problem of 
deri™g sCgtyl?ghtREE enriched magmas from light-REE depeleted source regions is common toterrestnal hot 
spots fuch as 8 Hawaii (8). Single-stage models require prohibitively small degrees of partial melting (<1 «). 
multi-stage melting models have been invoked to spread the REE fractionation over two or more steps (e.g. 9). 
Some son of multi-stage melting process thus seems necessary to explam the Nakhla parent magma composition 
wS the condition that garnet not have been a residual phase in the last stage of melting^ addition to low-Al the 
Nakhla parent magma also had an unusually high concentration of CaO (-14 wt% (2)). The combination of low 
mXgTcTXS either th.t the route. ». dominated by augite or dta. C0 2 . whteh ta, the potemtel o 
drastically increasing the CaO content of melts coexisiting with olivine and pyroxene, fluxed the melting 
pressures > 25 kb (2). Since partial melting of pyroxene dominated sources produces small negative Eu anomalies in 
die liquid (10) and since there is no evidence of such an anomaly in Fig. 2 A, the I^esence of C0 2 m mar 1 
melting processes must be seriously considered.C0 2 is also an effective earner of light REE (11), so C0 2 may have 

affected both trace and major elements during melting. __ * t vn* Tf , M 

Fig IB illustrates very different REE patterns for the shergottites. The Shergotty pattern is the 70% ICM 

model taken from (12). The EETA79001A pattern is the bulk rock analysis of (13). EETA79001A is a fine-grained 
basaltic rock with 10-15 % mafic xenocrysts (14). These xenocrysts will likely have only a minor dl ““ g 
incompatible elements, so the pattern in Fig. IB is believed to close to. albeit slightly lower and s^Pe^ than 
true parent liquid pattern. The crystallization ages of the shergottites are controversial because of variable sho k 
effects on the isotopic systems and consequentlythe values of ENd model dependent. The values shown m Fig. 

IB are consistent with the 180 m.y. age advocated by (15). This age is chosen here because only the younger ages, 

which yield £ Nd > 0. are petrologically reasonable, and because the 350 m.y. age reported by (16) has been shown 
to be a mixing line (12). Given these qualifications, the low-Al content, the depleted light-REE pattern, and +£ Nd of 
EETA79001A have a straightfoward explanation: partial melting of a ^w-^ source region with A^QOOlT^afma 
REE depletion. In this regard, the source region was similar to that of Nakhla although the EETA79001A magm 
genesis was apparently much simpler. The ENd values for Nakhla and EETA79001A are much higher than typical 
terrestrial basaltic values, but are more typical of lunar mare basalts. This similarity suggests that Mars was more 
like the Moon in its ability to maintain long term isotopic heterogeneities m its mantle. Lack of crustal recycling 
on Mars and/or less vigorous mantle convection than the Earth are probable explanations , Q(V ., . . 

Given the similarity of mineral compositions in Shergotty to those m the groundmass of EETA7900 , 
is likely that their parent magmas lay along similar liquid lines of descent, as suggested by Fig. 1, and hence they 
were derived from similar primary magmas and source regions. If so, then REE pattern and ENd of die Shergotty 


202 


MAGMATIC PROCESSES ON MARS: Longhi J. 


Shergotty parent magma in Fig. 2B are readily explicable as those of a magma derived from a depleted source 
region like EETA79001A, but subsequently contaminated by a low-temperature, long-term, light-REE enriched 
component. The slight U-shape in the light REE is especially indicative of such a contamination. This component 
probably is crustal, but whether it is older basalt, like the Nakhla parent magma (Fig. 2A), or 'granitic* is not 
clear; the physics of assimilation favors an evolved composition with a low melting point, however. One thing 
that is clear is the absence of a negative Eu-anomaly in the Shergotty REE pattern. Consequently, this crustal 
component was unlike lunar KREEP, which has a prominent negative Eu-anomaly (17). 

Fig. 2A contains REE concentrations for the Nakhla parent ('closed system') and the bulk data for 
EETA79001A taken from Fig. 1 plus additional calculated and measured concentrations of some high -field-strength 
elements (HFSE) arranged in order of incompatibility. Fig. 2A shows that there are complementary anomalies for 
Ta, Hf, and Zr in the Nakhla and EETA79001A patterns, thus supporting the hypothesis that the shergottites were 
generated by remelting a source depleted in a Nakhla-like component. Fig. 2B schematically illustrates typical 
incompatible element patterns for basalts from terrestrial oceanic islands (OIB), volcanic arcs (ARC), and mid- 
ocean ridges (MORB). Despite the fact that there is no evidence of plate tectonics on Mars, the Nakhla parent 
magma pattern appears more similar to the ARC pattern than to OIB or MORB. This similarity is probably due to 
similar fractionations of the REE from the HFSE during transport by a C02-rich vapor phase, rather than similar 
tectonic styles. 


REFERENCES: (1) McSween, H. Y. (1985) Revs. Geophys., 23, 391-416. (2) Longhi, J. and Pan V (1989) 
Proc. Lunar Planet. Sci. Conf. 19th, p. 451-464. (3)Toulmin, P„ HI, Baird. A. K.. Clark. B. C., Keil K Rose H 
L, Jr Christian, R. P„ Evans. P. H.. and Kelliker, W. C. (1977) J. Geophys. Res. 82, 4625-4634. (4) Bertka’ 
C. M. and Holloway. J. R. (1988) Proc. Lunar Planet. Sci. Conf. 18th, 723-739. (5) Treiman A H (1986) 
CMBWcWm - Acta - 50 - 1061-1070. (6) Nakamura, N.. Unruh, D. M.. Tatsumoto, M. and Hutchinson R 
(1982) Geochim. Cosmochim. Acta, 46, 1555-1573. (7) McKay G.. Wagstaff J.. and Yang S.-R. (1986) Geochim. 
Cosmochim. Acta, SO, 927-937 . (8) Chen C.-Y. and Frey F.A. (1985) J. Geophys. Res., 90, 8743-8768 (9) Ribe 
N.M. (1988) Earth Planet. Sci. Lett., 88, 37-46. (10) Shearer C.K. and Papike JJ. Proc. (1989) Lunar and Planet 
Set. Conf. 20th, m press. (11) Wendlendt R.F. and Harrison WJ. (1979) Contrib. Minerl. Petrol., 69, 409-419. 
(12) Lundberg L.L, Crozaz G„ McKay G.. and Zinner E. (1988) Geochim. Cosmochim. Acta, 52, 2147-2163. (13) 
Burghele A., Dreibus G„ Palme H., Rammensee W., Spettel B„ Weckwerth G„ and Wanke H. (1983) Lunar 

(14) McSween H Y - 311(1 Jarosewich E. (1983) Geochim. Cosmochim. Acta, 47, 
).„ 0 ':. 15 ' 3 ' (15 ) J ° nes J H - (1986) Geochim. Cosmochim. Acta, 50, 969-977. (16) Jagoutz, E. and Wanke H 
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Figure 1. REE in Nakhla 
and shergottite parent 
magma compositions. 
Nakhkla - calc.,this study; 
Shergotty - calc (12); 
EETA79001A - bulk (13). 




Fig. 2 A. Calculated 
REE and HFSE in o 

Nakhla parent • 

magma and bulk §“ 

EETA79001A. B. ■ 

Typical (schematic) 
patterns for 
terrestrial basalts. 
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VALLES MARINERIS, MARS: ARE PIT CHAINS FORMED BY EROSION AND TROUGHS 

BY TECTONISM?; B.K. Lucchitta, R.A. Balser, and L.M. Bertolini, 

U.S. Geological Survey, Flagstaff, AZ 86001. 

The origin of the Valles Marineris remains controversial. Erosional 
[1], tectonic [2], and hybrid processes [3] have been proposed. Erosional 
processes appeared attractive because a morphologic continuum was thought 
to exist from pit chains of probable erosional origin to larger troughs. 
Schultz [4], however, refuted the existence of this continuum. To clarify 
these contradictions, we compared the widths and depths of pit chains and 
troughs and found that the features do not form a continuum. Rather, 
results are consistent with the hypothesis that pit chains formed by 
surficial collapse and troughs by deeper seated and coherent failure. 

We classified by inspection all pit chains and linear depressions in 
the Valles Marineris region into six morphologic categories: (1) pit 

chains (linear arrays of small pits), (2) floored chains (arrays of pit 
chains having flat or hummocky floors), (3) scalloped troughs (wider linear 
depressions with scalloped wall segments), (4) narrow troughs (depressions 
of intermediate width with straight wall segments), (5) wide troughs 
(broad, linear depressions), and (6) chaotic troughs (more irregular 
depressions displaying some channel morphologies). We drew topographic 
profiles across the classified depressions at each degree of longitude 
between long 45° and 90°, on the basis of 1: 2,000,000-scale topographic 
maps of MC 18 NW [5] and 18 NE and SE (work in progress). For each 
profiled depression, we determined the erosional width between plateau 
margins, the depth from the surrounding plateau level to the deepest part, 
and the geologic unit (modified from Witbeck et al. [6]) exposed at the 

deepest point. . 

Depths and widths are shown in Figs. 1 and 2. Fig. 1 shows 
measurements of all troughs. Fig. 2 is the same but with omissions of 
three types: (1) chaotic troughs, which are genetically linked to outflow 

channels; (2) troughs east of long 61°, which are transitional to chaotic 
troughs; and (3) troughs whose maximum depth is not likely to reflect the 
approximate structural depth because of thick fill from landslides or 

interior deposits. , , 

Results of the study, as deduced from the figures, are as follows: 

Pit chains, floored chains, and scalloped troughs lie along a 
continuously ascending trend that suggests a common origin. 

Surficial erosional collapse into linear subsurface voids or tension cracks 
[3] is compatible with this observation. The limiting depths of about 4 km 
may be due to a discontinuity at that depth or to restricted size of the 

underlying voids. . 

2 . Narrow and wide troughs form a continuum that suggests their 
formation by similar processes. The straight scarps bordering these 
troughs suggest control by faulting rather than surficial collapse. 

3. A gap in data occurs at widths of 20-35 km, separating pit chains, 
floored chains, and scalloped troughs from narrow and wide troughs. Only 
two transitional points are located within this gap. The gap suggests an 
abrupt change in physical conditions or processes. Apparently, deep- 
seated, more coherent failure was activated for troughs wider than about 

35 km. , 

4. Most troughs bottom out at 8-9 km regardless of width, perhaps 

implying a controlling discontinuity or limit in the amount of extension. 
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A PHYSICAL MODEL OF THE IMPACTED MARTIAN CRUST: HYDROLOGIC 
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Photogeologic evidence suggests a close association of many outflow [1,2] 
intermediate [3], and sapping [4] channels with regions affected by uplift, fracturing, and 
geothermal heating [1, 5]. The magnitude and diverse morphology of the channeling are 
impressive by terrestrial standards, especially because there is little evidence that significant 
rainfall has occurred on Mars [6]. By the end of the heavy bombardment, outgassed water 
probably had been emplaced in the upper crust [7]. Nearly all the channels originated rom 
or are associated with ancient cratered terrain, some of which has been resurfaced by Java 
flows and eolian processes. Because high impact fluxes during bombardment probably 
fragmented most of the materials that resurfaced the impacted surface, we believe that a 
simple physical and structural model of the impacted, ancient cratered terrain can 
reasonably describe physical properties fundamental to interpreting subsequent channel 

formation, as well as other geologic activity. • „ t _ 

On this basis, we developed a two-layer model for the impacted crust, an ejecta 

(impact breccia) zone overlying a zone of fractured basement rock [see 8, 9] The proposed 
ejecta zone is 1 to 2 km thick and consists of well-mixed, very poorly sorted impact breccia 
that grades into the fractured zone. The fractured zone is composed of meter-size or larger 
blocks of basement rock, and it extends to the depth of hthostatic closure (10 km or more, 
depending on pore-water pressure [10]). The porosity and permeability of these zones 
depend on the size distribution and packing of the breccia and fractured blocks. From the 
results of measurements from terrestrial impact and explosion craters laboratory analyses 
and experimentally verified theoretical models for porosity and permeability, we determined 
that in such an impacted crust, clast size and degree of sorting increase with depth. These 
trends cause (1) porosity to decrease with depth throughout the affected crnst, (-j 
permeability to increase with depth to the top of the fractured zone and then to decrease 
as fracture width and frequency decrease with depth, and (3) material strength to increase 
with depth from cohesionless to that of solid bedrock. Extrapolating from typical size 
distributions and packing among ejecta and fractured blocks, we estimate the porosity o 
the ejecta zone as 10-20 percent and the porosity of the upper part of the fractured zone 
as less than a few percent; the permeability of the ejecta zone would be less than 0.0 
darcys, and the permeability of the upper part of the fractured zone would be about 1000 

We recognize two major limitations of our model: (1) few terrestrial impact and 
explosion craters provide empirical data for the analysis of the effects of single and multiple 
impacts on rock structure, and (2) the subsurface stratigraphy of Mars is poorly known. 
Thus some of our basic assumptions are speculative, allowing for different possibilities 
regarding the nature of the Martian crust. In one alternative view [10], the crust is more 
heterogeneous than we have portrayed; thus its hydraulic properties may be more varied 
and impact fractures in the basement rocks, if filled, would appreciably lower the bulk 

permeability.^ ^ problems, we feel that our two-layer model is generally valid, 
because it can account for a great diversity of Martian geologic phenomena. Photogeologic 
evidence suggests mechanical discontinuities in the Martian crust at 1- to 3-km depths, 
which in some places may form the contact between the ejecta and the fractured zones. 
Generally, lower layers of the ancient crust have progressively greater resistance to erosion. 
Within the deeply eroded Kasei Valles, for example, erosional discontinuities in the 
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channels are evident at depths of 1.0 and 2.6 km. The upper discontinuity has been 
interpreted as the ejecta/basement contact [11], the interface between ice-laden and dry or 
wet regolith [12], or the zone between pristine and cemented regolith [13]; the lower 
discontinuity is consistent with the depth to the base of sapping channels along the Valles 
Marineris and has been interpreted as the ejecta/basement contact [12]. Such 
discontinuities are consistent with observed graben widths and collapse pits proposed to 
originate by collapse of cohesionless material into tension fractures at depth in basement 
rocks [12, 14]. 

The Chryse channels may be explained by high pore pressures and fluidization of the 
base of the ejecta zone by the fractured zone, causing enormous retrogressive debris flows 
[15]: after removal of the overburden, floods from the fractured zone could have easily 
eroded ejecta-zone breccia. Late-stage debris flows, sapping channels, and landslides along 
the channel margins completed the morphology seen today. If debris flows played a central 
role in the formation of the large outflow channels, their mobilization would require much 
less water [16] than the flood model [7], The breccia of the ejecta zone is poorly sorted 
and rich in clay-size particles, analogous to the clast distributions of terrestrial debris flows 
[9]. An explanation for the great runout distances of debris flows over very low gradients, 
however, has not yet been offered. 

Other geologic features consistent with our model include (1) outflow channels 
produced by catastrophic floods erupted from joints or faults that tap the fractured zone 
and act as high-volume conduits [13,17]; (2) sapping channels where the permeability of the 
ejecta zone was sufficient to promote sapping of exposed, saturated ejecta; (3) chaotic 
terrain, possibly formed by liquefaction of breccia at depths of hundred of meters to more 
than a kilometer [18]; (4) complex channels such as Nirgal Vallis, which can be explained 
as a combination of outflow along a tension fracture that taps the fracture zone (lower part 
of channel) and sapping from the ejecta zone (upper part of channel); and (5) high-latitude 
debris aprons and channels with modest wall slopes (e.g., Auqakuh Vallis) composed of ice- 
rich impact breccia of low yield strength [19]. 

Our physical model for the ancient, impacted Martian crust, coupled with local 
volcano-tectonic histories, enables a more complete understanding of the formation of many 
common geologic features on Mars. Previously, catastrophic flood models have not 
addressed the dependence of permeability on clast and fracture distributions and the 
cohesion of the cratered terrain material, and debris-flow models have not addressed the 
physical mobility of the material. Furthermore, the ejecta/fractured zone stratigraphy 
explains many observations of Martian erosional and structural discontinuities. 
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